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Abstract 

We calculate the important next-to-leading-order (NLO) contributions to the B — ► KK* 
<^ \ decays from the vertex corrections, the quark loops, and the magnetic penguins in the per- 

turbative QCD (pQCD) factorization approach. The pQCD predictions for the CP-averaged 
branching ratios are Br( B + -> K + K*°) « 3.2 x 10~ 7 , Br{ B + -» K°i^*+) « 2.1 x 10~ 7 , 
Br(B°/B° -► ^°F° + K°iP°) « 8.5 x 10~ 7 , Sr( 5°/^° -► ^ + -KT*~ + K~ K*+) « 1.3 x 10~ 7 , 
which agree well with both the experimental upper limits and the predictions based on the 
QCD factorization approach. Furthermore, the CP-violating asymmetries of the considered 
decay modes are also evaluated. The NLO pQCD predictions for Acp{B + — ► K + K ) and 
A7p(5+ -> i^* + K°) are ^(if+^*°) « -6.9% and JVJi r p (K* + T?) « 6.5%. 
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I. INTRODUCTION 

It is well-known that the experimental measurements and theoretical studies of the 
two body charmless hadronic B meson decays play an important role in the precision 
test of the standard model (SM) and in searching for the new physics beyond the SM 
PJ. For these decays, the dominant theoretical error comes from the large uncertainty 
in evaluating the so-called hadronic matrix element, (MiM2\Oi\B), where M\ and M 2 
are light final state mesons. The perturbative QCD (pQCD) approach [2] is one of the 
most popular factorization approaches [3j, |4[ being used to calculate the hadronic matrix 
elements. 

When compared with the QCDF or SCET factorization approaches, the pQCD ap- 
proach has the following three special features: (a) since the kr factorization is employed 
here, the resultant Sudakov factor as well as the threshold resummation can enable us to 
regulate the end-point singularities effectively; (b) the form factors for B — > M transition 
can be calculated perturbatively, although some controversies still exist about this point; 
and (c) the annihilation diagrams are calculable and play an important role in producing 
CP violation. 

Up to now, almost all two-body charmless B/B s — > M L \M 2 d ecay s have been calculated 
by using the pQCD approach at the leading order [5|, |y, |7|, la, |9|, [HJ, llll, ll2( . Very recently, 



some next-to-leading (NLO) contributions to B — > Kir and several B — > PV decay modes 



[131 . |14| have been calculated, where the Wilson coefficients at NLO accuracy are used, and 
the contributions from the vertex corrections, the quark loops and the chromo-magnetic 
penguin operator 0% g have been taken into account. As generally expected, the inclusion 
of NLO contributions should improve the reliability of the pQCD predictions. 

In a previous paper [10], the authors calculated the branching ratios and CP violating 
asymmetries of the B°/B° -► K°K*° ,K°K*°, K+K*- , K~K* + , and B + -► K+K*° and 



K K* + decays by employing the pQCD approach at the leading order. Following the 
procedure of Ref. [l3j, we here would like to calculate the NLO contributions to the 
B — > K*K decays by employing the low energy effective Hamiltonian and the pQCD 
approach. 

The remainder of the paper is organized as follows. In Sec. II, we give a brief discussion 
about pQCD factorization approach. In Sec. Ill, we calculate analytically the relevant 
Feynman diagrams and present the various decay amplitudes for the studied decay modes 
in leading-order. In Sec. IV, the NLO contributions from the vertex corrections, the 
quark loops and the chromo-magnetic penguin amplitudes are evaluated. We show the 
numerical results for the branching ratios and CP asymmetries of B — > K*K decays in 
Sec. V. The summary and some discussions are included in the final section. 

II. THEORETICAL FRAMEWORK 

The pQCD factorization approach has been developed and applied in the non-leptonic 
B meson decays [2j] for some time. In this approach, the decay amplitude is separated 
into soft, hard, and harder dynamics characterized by different energy scales (t,m b , M w ). 



It is conceptually written as the convolution, 

A(B - M X M 2 ) ~ fd 4 hd 4 k 2 d% Tr ^(^^(^^(^^(^^(^fe.fes,*)] , (1) 

where fcj's are momenta of light quarks included in each meson, and Tr denotes the trace 
over Dirac and color indices. C(t) is the Wilson coefficient, which includes the harder 
dynamics at larger scale than Mb scale and describes the evolution of local 4-Fermi 
operators from mw down to t ~ O(vAMg) scale, where A = Mb — vn^. The function 
H(ki, k 2 , k 3 , t) describes the four quark operator and the spectator quark connected by a 

hard gluon whose q 2 is in the order of AMb, and includes the 0( y AMb) hard dynamics. 
Therefore, this hard part H can be perturbatively calculated. The function $^ is the 
wave function which describes hadronization of the quark and anti-quark into the meson 
Mi. While the function H depends on the processes considered, the wave function &Mi is 
independent of the specific processes. 

In the B meson rest-frame, it is convenient to use light-cone coordinate (p + ,p~ , p T ) to 
describe the meson's momenta, 

p ± = ^=(p°±p i ), and p T = {p\p 2 ). (2) 

Using these coordinates the B meson and the two final state meson momenta can be 
written as 

M B , s M B , , M B , , 

P b = -|(1,1,0t), P^ = -|(l,7i„0 T ), P i , = -|(0,l-4„OT), (3) 

respectively, here r#* = m^*jMB- The light meson (K) mass has been neglected. For the 
B — > K*K decays considered here, only the vector meson's longitudinal part contributes 
to the decays, and its polarization vector is e^ = Jff - (1, — r 2 K ,,^). Putting the anti- 
quark momenta in B, K* and K mesons as k\, k 2 , and k%, respectively, we can choose 

fci = (ziPi + ,0,ki T ), A; 2 = (x 2 P 2 + ,0,k 2T ), k 3 = (0,x 3 Pi,k 3 T)- (4) 

Then, the integration over k± , k 2 , and k 3 in eq.(pQ) will lead to 

A(B — > KK*) ~ dxidx 2 dx 3 bidbib 2 db 2 b 3 db 3 

■Tr [C(t)$ B (x 1 ,6 1 )$^(a; 2 ,6 2 )$^(x3,63)^(x i ,6 i ,t) > S t (a; i )e- 5(i) ] , (5) 

where foj is the conjugate space coordinate of k^, and t is the largest energy scale in func- 
tion H(xi,bi,t). The large logarithms (Inmyy/t) coming from QCD radiative corrections 
to four quark operators are included in the Wilson coefficients C(t). The large double 
logarithms (In 2 xi) on the longitudinal direction are summed by the threshold resumma- 
tion 15j], and they lead to St(xi) which smears the end-point singularities on X{. The last 



term, e~ s<yt \ is the Sudakov form factor which suppresses the soft dynamics effectively [2J, 
Thus it makes the perturbative calculation of the hard part H applicable at intermediate 
scale, i.e., M B scale. 



A. Wilson Coefficients 

For B — > KK* decays, the related weak effective Hamiltonian H e ff with b —>■ s tran- 
sition can be written as |l6j 

n eff = % { E V #K [ c M° q M + c^cM - v tb v; d f; c^) o^) \ . (6) 

* \q=u, c i=3 J 

with GV = 1.16639 x 10~ 5 GeV~ 2 is the Fermi constant, and Vij is the CKM matrix 
element, Ci(fx) are the Wilson coefficients evaluated at the renormalization scale // and 
Oi are the four-fermion operators. For the case of b — » d transition, simply make a 
replacement of s by d in Eq. (JHJ) and in the expressions of Oi operators, which can be 



found easily for example in Refs. [lOj, [16 



In PQCD approach, the energy scale "£" is chosen at the maximum value of various 
subprocess scales to suppress the higher order corrections, which may be larger or smaller 
than the mj, scale. In the range of t < rrib or t > mb, the number of active quarks is Nf = 4 
or Nf = 5, respectively. For the Wilson coefficients Ci{[i) and their renormalization group 
(RG) running, they are known at NLO level currently [16] . The explicit expressions of 
the LO and NLO Ci(mw) can be found easily, for example, in Refs. p, ll6j . 



When the pQCD approach at leading-order are employed, the leading order Wilson 
coefficients Ci(mw), the leading order RG evolution matrix U(t, m)^ from the high scale 
m down to t < m ( for details see Eq. (3.94) in Ref. [16J), and the leading order a s (t) are 
used: 

47T 
° s ^ = o 1 U2/A2 T> ( ? ) 

A) In [t 2 /A 2 QC D\ 

where (3 = (33 - 2JV»/3, K { § CD = $.22hGeV and h!§ CD = 0.287 GeV. 

When the NLO contributions are taken into account, however, the NLO Wilson coef- 
ficients Ci(mw), the NLO RG evolution matrix U(t,m,a) (for details see Eq. (7.22) in 



Ref. [16]) and the a s (t) at two- loop level are used: 



„ m _ 4vr f A ln[ln[t 2 /A^ CD ]] 

where /3 = (33 - 2JV»/3, /3i = (306 - 38JV»/3, A^ D = 0.225 GeV and A^ D = 0.326 
GeV. 

By using the input parameters as given in the Appendix, it is easy to find the numerical 
values of the LO and NLO Wilson coefficients Cj(m&) for m^ = 4.8 GeV, as listed in Table 

m 



B. Wave Functions 

The B meson is treated as a heavy-light system. For the B meson wave function, since 



the contribution of (fi B is numerically small [17] , we here only consider the contribution 



TABLE I: The numerical values of the LO and NLO Wilson coefficients Ci(m b ), C7 7 (m;,) and 
C Sg (m b ). 



Ci(m b ) 


C x 


c 2 


c 3 


c± 


c 5 


c 6 


LO 


-0.2812 


1.1246 


0.0130 


-0.0278 


0.0080 


-0.0343 


NLO 


-0.1747 


1.0774 


0.0125 


-0.0330 


0.0094 


-0.0393 


Ci(m b ) 


C 7 /a 


C s /a 


C 9 /a 


Cio/a 


C77 


C 8 g 


LO 


0.1338 


0.0514 


-1.1459 


0.2865 


-0.3109 


-0.1481 


NLO 


-0.0032 


0.0305 


-1.2760 


0.2553 


-0.3016 


— 



of Lorentz structure 
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V2Nc 



i^B + m B )j 5 4> B (k 1 ) 



with 



(9) 



^(x, b) = N B x 2 (l — x) 2 exp 



Ml x 2 1 



M 



M) 1 



(10) 



where u b is a free parameter and we take u b = 0.4 ± 0.04 GeV in numerical calculations, 
and N B = 101.445 is the normalization factor for u b = 0.4. 

The K and K* mesons are all treated as a light-light system. The wave function of K 
meson is defined as 18 



*k{P,x,0 



1 



V2N ( 



=75 Wk( x ) + ™q<I>k(x) + (m, 



— v ■ n) 



£ 



x 



'ir 



c 



where P and x are the momentum and the momentum fraction of K, respectively. The 
parameter ( is either +1 or —1 depending on the assignment of the momentum fraction 
x. For the considered B — > KK* decays, K* meson is longitudinally polarized, and only 
the longitudinal component 0^* of the wave function contribute |18| 



& 



1 



V2Nc 



{h [m K ,(j) K ,(x) + ^k*</>k*(x)] + m K *4>K*(x)} 



(12) 



where the first term is the leading twist wave function (twist-2), while the second and 
third term are sub-leading twist (twist-3) wave functions. 

The expressions of the relevant distributions functions are the following 18J: 



<j> A K {x) 
<f> p K {x) 
f K {x) 



I 



K 



Ik 



:Qx(l — X) 



Krfi/2 



K r ,3/2, 



K'^3/2, 



1 + of Cf\t) + a«C*'\t) + af Cr(t) 



^1/2, 



2y/2N c 
2y/2N c 



l + (30 m --p 2 K )C^(t)-3 



9 
W4j + ^P* (1 + 6a 2 ) 



yl/2 



cr® 



-2 ^ o i, : 



I + 6(5773 - -r/3^3 - ^Pk " f /&a?)(l " 10x + 10x 2 ) 



(13) 

,(14) 
(15) 



with the mass ratio px = tuk/itiok, and T]3 = 0.015, u = —3.0. Since the uncertainties 



of the currently available Gegenbauer moments [19| are still large, we vary the value of 
a\ and a% by 100%, i.e. af = 0.17 ± 0.17, af = 0.115 ± 0.115, but keep a% = -0.015, 
because the theoretical predictions are insensitive to o^ K . 



as: 



The twist-2 DAs for longitudinally polarized vector meson K* can be parameterized 

<f> K .(x) = ^=6^(1 - x) [l + a XK *Cl /2 {t) + a 2K ,Cl'\t)\ , (16) 

where fx* = 200 MeV is the decay constant of the vector meson with longitudinal polar- 
ization, and the Gegenbauer moments are cl\k* — 0.03 ± 0.03, a 2 K* — 0.11 ± 0.11. As for 
the twist-3 DAs <fi s K * and </4:*, there is no recent update associated with those updates for 
twist-2 DAs, we adopt their asymptotic form: 



>f K ,( x ) = M=(l-2x), ^(x) = l%=(2x - \ Y (17) 



2y/2N~ c "~ ~~ n ™~ v ~' 2^W C [/J 



At last the Gegenbauer polynomials C^(t) are given as: 

Cl /2 {t) = i(3t 2 - 1), C\ /2 {t) = 1(3 - 30t 2 + 35t 4 ) 
Cl /2 {t) = 3t, C 2 3/2 (t) = ^(5t 2 -l), 
C 4 3/2 (t) = ^(l-l4t 2 + 21t 4 ), 

o 

with t — 2x — 1. 



III. DECAY AMPLITUDES AT LEADING ORDER IN PQCD APPROACH 



The B — > KK* decays have been studied previously in Ref. [10J by using the lead- 
ing order pQCD approach. In this paper, we focus on the calculations of some NLO 
contributions to these decays in the pQCD factorization approach. For the sake of com- 
pleteness, however, we firstly recalculate and present the relevant LO decay amplitudes 
in this section. 

At the leading order, the relevant Feynman diagrams for B° — > K*°K , K°K , B° — > 
K + K*-,K-K* + , and B + -»■ K + K*°, K* + K° decays have been shown in Figs. d [2] and O 

As illustrated by Fig. [TJ both B° and B can decay into K*°K and K°K simulta- 
neously. Besides of the eight Feynman diagrams in Fig. [TJ other four Feynman diagrams 
can be obtained by connecting the left-hand end of the gluon line to the lower d quark 
line inside the B° meson for (e) and (f), or to the lower s or d quark line for (g) and (h). 

For B° — > K K and K°K decays, only the operators O3-10 contribute via penguin 
topology. Its is a pure penguin mode with only one kind of CKM element £j = V t * b V t d, 
and therefore there is no CP violation for such decays at leading order. 

For B° — > K*°K , K°K decays, we firstly consider the case of B — > K*° transition 
where K*° meson takes the spectator d quark. For the (V — A)(V — A) operators, the 
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FIG. 1: Diagrams contributing to the B° — > K*°K decay. From diagram (a) and (b), the form 
factor Aq^ k or F B ^ K can be extracted. Other four Feynman diagrams can be obtained by 
connecting the left-hand end of the gluon line to lower d quark line inside the B° meson for (e) 
and (f), while to the lower s or d quark line for (g) and (h). 

decay amplitude corresponding to Figs. ([T])a and ([T])b can be written as 

FeK* — 4\/2G ' Fl\C 'f J "k^B \ dX\dX2 / bidbib 2 db 2 (j)B(Xl,bi) 

Jo Jo 

x { [(1 + x 2 )(J)k* (x 2 ) - (1 - 2x 2 )r K * (cj) s K * (x 2 ) - 4>k* (x 2 ))] 

■E e {t a ) h e (x 1 ,x 2 ,b 1 ,b 2 ) -2r K *<\) s K *{x 2 ) ■ E e (t' a ) h e (x 2 ,x 1 ,b 2 ,b 1 )} , (19) 

where r K * = m K */mB, Cp = 4/3 is a color factor. The evolution factors E e (t a ) and the 
hard functions h e are displayed in Appendix lAl 

For the (V - A)(V + A) and (S - P){S + P) operators, we find 



F P1 - - F is* 

r eK* ~ r eK* 



(20) 



/"l /"OO 

FeK* = 8v2GFnCFfK'm' B / dx\dx 2 \ bidbib 2 db 2 (j>B(xi) 

Jo Jo 



x {-r K [4> K *(x 2 ) - r K *{{2 + x 2 )<j) a K *{x 2 ) + x 2 <t> t K *{x 2 ))\ 

■E e (t a ) h e (x 1 ,x 2 ,b 1 ,b 2 ) + 2r K *r K (t) s K ,{x 2 ) ■ E e (t' a ) h e (x 2 ,x 1 ,b 2 ,b 1 )} . (21) 



For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 



involved. The decay amplitudes are 

M e K* = -^Gf^Cfvo^ / dx\dx 2 dx 3 \ bidbib 3 db 3 (f>B(xi,bi)(j) K (x 3 ) 
v3 Jo Jo 

x { [r K *x 2 (0x*fe) + 0k-.(x 2 )) + (1 - x 3 )(f) K *(x 2 )] 

•E' e (tb)h n (x 1 ,x 2 , 1 - x 3 , 6i, 6 3 ) + E' e (t' b )h n (x 1 ,x 2 , x 3 , b ± , b 3 ) 

■ [-(x 2 + x 3 )(j) K *(x 2 ) - r K *x 2 (<j) s K *(x 2 ) - <f>K*(x 3 ))] } , (22) 

M-zk* — -j=GFT\CFvr^ B r K \ dx\dx 2 dx 3 \ b\db\b 3 db 3 4>b{x\, 61) x {[(1 — x 3 )4>k*{x 2 ) 

v3 Jo Jo 

■ Wfe) - 0k-(x 3 )) -r K *(l- x 3 ) ((j) s K * (x 2 ) + 4> l K * (x 2 )) ((J>k(x 3 ) - 0^(x 3 )) 
-r K *x 2 ((j) s K *(x 2 ) -$r*(s 2 )) (0^(x 3 ) + 0x(^3))] E' e (t b )h n (x 1 ,x 2 , l-x 3 ,h ,b 3 ) 

- [x 3 ■ (f)K*{x 2 ) (<Pk( x 3) + 4>k{^)) - r K *X 3 {(t>K*{x 2 ) + (J)k*(x 2 )) (0£(x 3 ) 

+4> T k{x 3 )) - r K ,x 2 {<p s K *{x 2 ) - ^ K *{x 2 )) (0xfe) - 4>k(x 3 ))] 
xE' e (t' b )h n (x 1 ,x 2 , x 3 , 61, 63)} , (23) 

Mg. = 0. (24) 

For the non-factorizable annihilation diagrams (e) and (f), again all three wave func- 
tions are involved. The decay amplitudes are 

M a K* = —j^G pT^C Ffn B / dx\dx 2 dx 3 I bidbib 3 db 3 <f)B(xi,bi) 

v3 Jo Jo 

x { [(1 - x 2 )0^*(x 2 )0^(x 3 ) + r K *r K {\ - x 2 ) (4> s K *{x 2 ) + 0^,(x 2 )) {<p P K (x 3 )- 
4> T k{x 3 )) + r K *r K x 3 (<f) s K *(x 2 ) - 0^*(x 2 )) (<££(x 3 ) + 0^(^s))] 
x^(i c )/i na (xi,x 2 ,x 3 ,6i,6 3 ) - E' a {i!^H na {xx-,x 2 , x 3 ,h,b 3 ) [x 3 4)K*{x 2 )(t>K{x 3 ) 
+4r K *r K (f) s K *(x 2 )<f) P : (x 3 ) - r K *r K (l - x 3 ) (4> s K *(x 2 ) + <^*(x 2 )) (0£(z 3 ) 
-^(^3)) - r x *r x x 2 (<j> 8 K .(x 2 ) - ^.(x 2 )) (^(x 3 ) + 0xfe))] } , (25) 



M" = ^ 



16 Z" 1 Z" 00 

—=Gf^Cfiti b / dx\dx 2 dx 3 / bidbib 3 db 3 4>B{xi,bi) 
v3 Jo Jo 

{ [(1 - z 2 )r*r*</4(z 3 ) (^(^2) + 0^-,(^ 2 )) - r K x 3 (t)K*{x 2 ) (<Pk(x 3 ) - <Pk( x 3))] 
xE' a (t c )h na {xi ) x 2) x 3 MM) - [~(x 2 + l)r K *(fr^(x 3 ) (4> s K *(x 2 ) +$r*(x 2 )) 
-r K (x 3 - 2)0 x *(x 2 ) (0£(s 3 ) - <t> T K {x 3 ))] E' a {t' c )ti na {x u x 2 ,x 3 ,b u b 3 )} ■ (26) 

M^* = —=GF^CFfn B I dx\dx 2 dx 3 \ b\db^ 3 db 3 $B{x\,b\)\\{x 2 — \) 
v 3 Jo Jo 

x0^(S 2 )0xfe) - 4r x r x *0^*(x 2 )0^(x 3 ) + r K r K *x 2 (4> s K *(x 2 ) + <^*(z 2 )) 

• (^(^s) - 4>k( x s)) + r K *r K (l - x 3 ) (<f) s K *(x 2 ) - <^*(z 2 )) (^(xs) + </>kM)] 

■E' a (t e )h na (x 1 ,X 2 ,X 3 ,b 1 ,b 3 ) + [ x 3<pK* {x 2 )4>k{x 3 ) + X 3 T K *r K (<j) S K »(x 2 ) +0* x *(x 2 )) 

• (0xfe) - 0xfe)) + r^*r^(l - x 2 ) (</> 8 K .(x 2 ) - <&.(z 2 )) (^£(x 3 ) + $r(x 3 ))] 
x^(OC(x 1>a ; 2 ,X3,6i,63)}- (27) 



The factorizable annihilation diagrams (g) and (h) involve only K* and K wave func- 
tions. There are also three kinds of decay amplitudes for these diagrams. F a x* is for 
(V-A)(V-A) 

rl poo 

F aK * = F[x„ =4V2G F 7rC F f B m i B / dx 2 dx 3 / b 2 db 2 b 3 db 3 {- [(1 - x 2 )0ic*(x 2 ) 

Jo Jo 

■<pi{x 3 ) + Ar K r K *(t) s K ,(x 2 )(p p K {x 3 ) ~ 2r K *r K x 2 (t) P K {x 3 ) ((j) s K *{x 2 ) + <&.(s 2 ))] 
•-B a (td)/ia(^3,l-^2,&3,&2) + [a;30x*fe)0^(^3) + 2r K r K »<\) s K , (x 2 ) 
■ (^fe) + 0^(^s)) + 2r K r K *x 3 (f) s K ,(x 2 ) (0£(z 3 ) - ^jfe))] 
x J B a (t , d )/ ia (l-X2,X3,6 2 ,M} , (28) 

^5* = -SV2G F 7rC F f B m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x { [r A '*(l - x 2 ) (0^.(x 2 ) - <t>K*{x2)) <p£{x 3 ) + 2r K <f) K *{x 2 )<j) P : {x 3 )] 
xE a (t d )h a (x 3 , 1 - x 2 , b 3 , b 2 ) 

+ [2r K *(j) s K «(x 2 )(j)K(x 3 ) +x 3 r K (j) K *(x 2 )((j) P : (x 3 ) + <p T K (x 3 ))] 
xE a (t' d )h a {l-x 2 ,x 3 ,b 2 ,b 3 )} . (29) 

For the case of B° —>■ K° transition where K° meson takes up the spectator d quark, 
as shown in Fig. [Tj, it is straightforward to find the decay amplitudes by following the 
same procedure as the case of B° —>■ K*° transition. 

pi poo 

F eK = lQirC F m B I dx\dx 2 \ b\db\b 2 db 2 $B{x\,b\) 
Jo Jo 

x { [(1 + x 2 )4>i(x 2 ) + (1 - 2x 2 )r K ((j) P K (x 2 ) - 4> T K (x 2 ))\ E e (t a )h e (x 1 , x 2 , h, b 2 ) 
+2r K (fr P : (x 2 )E e (t' a )h e (x 2 ,x 1 ,b 2 ,b 1 )} , (30) 

*£ = F eK} F p £ = 0. (31) 

16 f 1 f°° 

M e x = —j=Gp^C F m B \ dx\dx 2 dx 3 \ bidbib 3 db 3 (pB(xi,bi)(f>K*(x 3 ) 
v3 Jo Jo 

x { [-r K x 2 {<p P K (x 2 ) + 0£(x 2 )) + (1 - x 3 )0%(x 2 )] 
■E' e (t b )h n (x 1 , x 2 , 1 - x 3 , b x , 63) + E' e {t' h )h n {xi, x 2 , x 3 , h, b 3 ) 

■ [-(x 2 + x 3 )(J)k(x 2 ) + r K x 2 ((I>k(x 2 ) - (t) T K {x 3 ))] } , (32) 

16 f 1 f°° 

M[k = -=.G p'kC pva B r k* \ dx\dx 2 dx 3 \ b\db\b 3 db 3 (J>b{xx, &i) x { [(1 — x 3 )(p^(x 2 ) 

v3 Jo Jo 

■ (<t>K*{x*) - <f>K*(x 3 )) + r K {l - x 3 ) (4>k(x 2 ) + 4> T K {.x 2 )) (</> 9 K *(x 3 ) - <t>K*{xz)) 
+r K x 2 (4>k(x 2 ) - 4> T K {x 2 )) {<p s K * (^3) + <Pk* (^3))] E' e (t b )h n (xi, x 2 , l-x 3 ,h, b 3 ) 
- [x 3 ■ (f>i{x 2 ) {(j) s K *{x 3 ) + 4>k*{x 3 )) + r K x 3 (<j) P K {x 2 ) + 4> T K {x 2 )) {(j) s K *{x 3 ) 
+<f>K*{xa)) +r K x 2 (4>k(x 2 ) - <f K {x 2 )) (4> s K *(x 3 ) - (Pk*(x 3 ))} 
xE' e (t' b )h n (x 1} x 2} x 3 MM)} ■ (33) 



1 f\ /*! /*oo 

M a K = —j=G fttC ' Fin's / dx\dx 2 dx 3 / bidbib 3 db 3 4>B{xi-,bi) 
v3 Jo Jo 

x { [(1 - x 2 )(J)k(^2)<Pk*(x 3 ) - r K r K ,(l - x 2 ) (V>£(x 2 ) + 4>k(x 2 )) ((p s K ,(x 3 )- 
0x*(%)) - r Kr K *x 3 (0xfe) - 4> T k{x 2 )) (<f) s K *(x 3 ) + 0^.(x 3 ))] 

xE' a (t c )h na (xi ) X 2) X 3 , 61, 63) - E' a (t' c )h' na (x 1 ,X 2 , X 3 , 61, 63) [^3</4(^2)0K*(^3) 

-4r x r^.0^(x 2 )0x*(^3) - r K r K *(l - x 3 ) (#£(x 2 ) + ^(x 2 )) (0x*fe) 
-0a-.(x 3 )) + ^r^xa (^(x 2 ) - 0xfe)) (0x-.(x 3 ) + <&.(* 3 ))] } , (34) 

M^ = —^Gp^CFfn^B \ dx\dx 2 dx 3 I bidbib 3 db 3 <f>B(xi,bi) 
v3 Jo io 

{ [(1 - x 2 )r*r*0*r*(x 3 ) (#£(x 2 ) + ^(^2)) - r^*x 3 ^(x 2 ) (^.(x 3 ) - 0a-.(x 3 ))] 
xE' a (t c )h na (xi,x 2 , x 3 ,bi,b 3 ) - [(x 2 + l)r K <p K *(x 3 ) (0£(z 2 ) + </>#■ (^2)) 
-r^,(x 3 -2)0^(x 2 ) (0^*(x 3 ) -0a-.(x 3 ))] ^(^)/i(, (xi,X2,X3,6i,6 3 )} , (35) 

M^ = — -=G fkC Fifn B / dx\dx 2 dx 3 \ b±dbib 3 db 3 4>b(xi : 61) {[(#2 — 1) 
v 3 Jo io 

x0^(S 2 )0x*(^3) + 4r^.r x 0^(S 2 ) 0^.(^3) - r K *r K x 2 ((j) P K {x 2 ) + 0^(x 2 )) 

• (0x*fe) - <Pk*(x3)) ~ r K r K *{l - x 3 ) {<p P K {x 2 ) - <f K {x 2 )) (<p s K *(x 3 ) + ^(x 3 ))] 
•^(i e )/i n a(^i,^2,a; 3 ,fci,fe 3 ) + [a:30^fe)0x*(%) - x 3 r K r K * (0£(x 2 ) + 0xfe)) 

• (0x*fe) - 0x*fe)) - r x r^.(l - x 2 ) (^(x 2 ) - ^(x 2 )) {^k*(x 3 ) + <&.(s 3 ))] 
xE' a (t' e )h' na (x u x 2 ,x 3 A,h)}- (36) 



/*1 /»00 

-^.^ = F^ = A\Z2G F ^C F fBin B / dx 2 dx 3 I b 2 db 2 b 3 db 3 {- [(1 ~ x 2 )(f)^(x 2 ) 

Jo Jo 

•<t>K*{x 3 ) - 4r^*r^0^(x 2 )0^-.(S3) + 2r K r K ,x 2 (p s K *(x 3 ) (0^(x 2 ) + ^(x 2 ))] 
•E a (t d )h a (x 3 , 1 - x 2 ,b 3 ,b 2 ) + [x 3 0^(x 2 )0^»(x 3 ) - 2r K *r K (t>K{x 2 ) 
• (<^*(x 3 ) + 4>K'(x 3 )) ~ 1r K ,r K x 3 (\) P K {x 2 ) (<f) s K *(x 3 ) - <&.(s 3 ))] 
x£' (i' d )/i a (l -z 2 , x 3 , 62,63)} , (37) 

F^ = 8v2G p^C p fs m B / dx 2 dx 3 / b 2 db 2 b 3 db 3 

Jo Jo 

x { [r K (l - x 2 ) {<p P K (x 2 ) - 4>k(x 2 )) <Pk*(x 3 ) - 2r K *<f)£(x 2 )<f) s K *(x 3 )] 
xE a (t d )h a (x 3 , 1 - x 2 , b 3 , b 2 ) 

+ [^■ k <Pk{.x 2 )(J) K *(x 3 ) - x 3 r K *(j)^(x 2 ) {4> s K *{x 3 ) + 0' X *(S 3 ))] 
xE a (t' d )h a (l-x 2 ,x 3 ,b 2 ,b 3 )} . (38) 
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FIG. 2: Diagrams contributing to the -B — > K + K* {K* + K ) decays. 



Combining the contributions from different diagrams in Fig. [TJ the total decay ampli- 



*n 



tude for B — > K K and K K decay can be written as 



M(B° - K*°T? + K°K*°) = -& \ (F eK + F; K ) ( a 4 - Ja 10 ) + F e % [ a, - -a s 



1 



+ 



(M eK + M* eK ) fc 3 - l -C^j + (M* + MS.) (<7 5 - ic 7 ) 



+ (M a * + M aK .) IC 3 + 2C 4 - -C 9 - C 10 

+ (M* + M&) (a - ^r) + (Af£ + M^,) (2C 6 - C s 



+F aK I 2a 3 + a 4 + 2a 5 - a 7 - ag - -a W/ 

+F aK * (2a 3 + a,-a 9 - ^o 10 )J + (F£ + F%,) (a 6 - ^a 8 J | 



(39) 



For B° — > K + K* (K* + K ) decays as illustrated in FigfSJ only annihilation diagrams 



-ftO 



contribute at leading order. Again, both B and B mesons can decay into the final state 
K + K*~ and its charge-conjugate state K~K* + . For B + — ► K+lT and K* + T( decays, 



II 
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K+(K* 




K+(K* 



s 



u K+(K 



*+\ 
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(d) 




-s 
-u 





u tf+(tf*+) 




FIG. 3: Diagrams contributing to the B + — ► K + K (K* + K ) decays. 



as shown in Fig. [3], the factorizable emission diagram, the non-factorizable spectator and 
annihilation diagrams contribute simultaneously. 

Following the same procedure as for B° — > K*°K /K°K decays, we find the total 
decay amplitude for the later two decay modes: 

M (B° -► K + K*- + K~K* + ) = £ u [(M aK + M aK *\ C 2 + (F aK + F aK «)a 2 ] 

-6 l(M aK + M aK «) Uc, + ho x ^j + (M% + M%) Uc, + ~C 8 ) 



1 l 

+ (F aK + F aK *) ( 2a 3 + 2a 5 + -o 7 + ~a 9 



(40) 



jV((5- 



ir+ir 



■o*. 



£« (M a ^ Ci + F aK ax)-£t{ FeK ( «4 - -aio 



r/'i in , ^ \ , 7^ /_ , _ \ , t^P 2 



+M eX ( C 3 - -C 9 + M^ Cs - ^C 7 ) + M ai , (C 3 + C 9J 



+M* (C 5 + C 7 ) + F ai ^ (a 4 + aio) + i 7 ^ K + a 8 ) } 



(41) 
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M{B + -► K*+IC) = e„ {M aK * d + F ai ,* oi) 



-6 <{ i 7 ^* ( 04 - g °10 ) + F eK* ( a 6 ~ ^ a 8 



+M eiC * ( C 3 - ^C 9 ) + M e p ^ ( C 5 - ^-C 7 



+M aK , (C 3 + C 9 ) + M a p ^ (C 5 + C 7 ) 

+i^* (a 4 + aio) + i^. (a 6 + a 8 ) } , (42) 

In the decay amplitudes of Eqs. (|3"9l - (l4"2"j) . the coefficients a,, the standard combination 
of the Wilson coefficients Ci, have been defined as usual 

n . C\ _ C 2 

d\ — 02 + — , 02 — Oi + — , 

C- 

a>i = Ci H — , for 2 = 3,5,7,9, 

o 

. = Q + ^zi 5 for j = 4,6,8,10. (43) 

Based on these decay amplitudes, the leading order pQCD predictions for the branching 



ratios and CP violating asymmetries of the considered decays can be calculated [10 



IV. NLO CONTRIBUTIONS TO B ^ KK* DECAYS IN PQCD 

The power counting in the pQCD factorization approach [13j is different from that in 
the QCD factorization [3J. When compared with the previous LO calculations in pQCD 
lOj, the following NLO contributions will be included: 



1. The LO Wilson coefficients Cj(mw) will be replaced by those at NLO level in NDR 
scheme [16J]. As mentioned in last section, the strong coupling constant a s (t) at 
two-loop level as given in Eq. flSJ), and the NLO RG evolution matrix U(t, m, a), as 



defined in Rcf. 161, will be used here: 



0. 
UOmi, m 2 , a) = U(m 1 ,m 2 ) + —R(m 1 ,m 2 ) (44) 

47T 

where the function U(mi, m 2 ) and R(mi, m 2 ) represent the QCD and QED evolution 



and have been defined in Eq. (6.24) and (7.22) in Ref. [16j. We also introduce a 



cut-off A cu t = 1 GeV for low energy scale in the final integration. 

2. The NLO contributions to the hard kernel H, including the vertex corrections, the 
quark loops, and the magnetic penguin |l3 



A. Vertex corrections 



The vertex corrections to the factorizable emission diagrams, as illustrated by Fig. (J2J), 
have been calculated years ago in the QCD factorization appeoach[3j, |20j. According to 
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B 



FIG. 4: NLO vertex corrections to the factorizable amplitudes. 



Ref . [13] , the difference of the calculations induced by considering or not considering the 
parton transverse momentum is rather small, say less than 10%, and therefore can be 
neglected. Consequently, one can use the vertex corrections as given in Ref. [20] directly. 
The vertex corrections can be absorbed into the re-definition of the Wilson coefficients 
aj(/i) by adding a vertex-function Vi(M) to them |3L l20l 



(kip) 



«*M + ^CV^#VKM), for i 



a,0) + 



47T 



1,2; 



4tt 



-C 



N r 



Vj(M), for j 



10, 



(45) 



where M is the meson emitted from the weak vertex. When M is a pseudo-scalar meson, 
the vertex functions Vi(M) are given ( in the NDR scheme) in Refs. [13j, |20j: 



Vi(M) 




fi dx(ftff(x)g(x), 
f^dx(^(x)g(l-x), 



for 



J Q dx(f>M(x)h(x) 



i = 1-4,9,10, 
for i = 5,7, (46) 

for i — 6, 8, 



where fu is the decay constant of the meson M; $%{x) and 4>m( x ) are the twist-2 and 
twist-3 distribution amplitude of the meson M, respectively. For a vector meson M, 



wwm) ^ s re pl ace d by (J)m{4> s m) and fu by f'lf in the third line of the above formulas. 



KiiKi) is replaced by 0m(0m) and hi by fu 

The hard-scattering functions g(x) and h(x) in Eq. PBT) are: 

l-2z 



^ia;j 



In a; — m 



x 



+ 



2Li 2 (x) — In x + 



21nx 



.r 



(3 + 22vr)lnx 



x 



/i(x) = 2L%2{x) — In x — (1 + 2«7r) lnx — [x ^ 1 — x) 



x) . (47) 

(48) 

where Li 2 (x) is the dilogarithm function. As shown in Ref. [l3j, the /^-dependence of 
the Wilson coefficients aj(/i) will be improved generally by the inclusion of the vertex 
corrections. 



B. Quark loops 

The contribution from the so-called "quark-loops" is a kind of penguin correction with 
the four quark operators insertion, as illustrated by Fig. (J5J). In fact this is generally 
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(a) 
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n-F^O , T ^n-?7*0 



--ftO 



FIG. 5: Quark-loop diagrams contributing to B° -> K*°K + K°K , B+ -> iT*+^ and 
K + K decays. 



called BSS mechanism [21[], which plays a very important role in CP violation. We here 
include quark-loop amplitude from the operators Oi j2 and 3 _ 6 only. The quark loops 
from O7-10 will be neglected due to their smallness. 

For the b — > d transition, the contributions from the various quark loops are described 
by the effective Hamiltonian H^K 



H (4) 
U eff 



E E % Wist ° {q) ^ /2) ^ (1 ~ 7s) T%) (?VT V) ' (49) 



q=u,c,t ql 



where I 2 being the invariant mass of the gluon, which connects the quark loops with the 
q'q pair as shown in Fig. [5j The functions C^(/x, I 2 ) can be written as 



C^\^f 



G {q \fi,l 2 



C2M, 



(50) 



for q = u,c and 



G (s \fx,f 



CM+ J2 G«'\»,l 2 )[C 4 (v) + CM}. (51) 

qf/=u,d,s,c 



The integration function G^ q '(fjL, I 2 ) for the loop of the quarks q = (u, d, s, c) is defined as 
13 



G {q \^l 2 ) = -4 I dxx{l-x)\n 
Jo 



m 2 — x(l — x)l 2 

V 2 



(52) 



where m q is the quark mass. The explicit expressions of the function G*(jU, I 2 ) after the 
integration can be found, for example, in Ref. [13] . 

It is straightforward to calculate the decay amplitude for Fig.(j5])a and (jHDb. For the 
case of B — > K* or B — > K transition, we find the corresponding decay amplitude M^l K 
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and Mjl K , with q — u,c, t, respectively; 

A pi PtX) 

Mk*k = — -j=GfC 2 f vti b I dx\dx 2 dx 3 I b\db\b 2 db 2 (f)B{xi,bi) 

v3 Jo Jo 

-r K « (1 - 2x 2 ) [<j> 8 K *{x 2 ) - 4>k*(x 2 )] <J>k{x 3 ) - 2r K <\> K «{x 2 )(f K {x 3 ) 
+2r K *r K [(2 + x 2 ) (f> s K *(x 2 ) + £ 2 <^,(x 2 )] (j) P K {x 3 )} 

■E^\t q ,l 2 )h e {x 2 ,x u b 2 M) 
+ {-2r K *(j) 8 K ,{x 2 )(j)i{x 3 ) + ±rK*r K 4>K*{x 2 )<j&{x z )} 

■E^(t' q J' 2 )h e (x u x 2 ,b u b 2 )} (53) 

and 

M-£ K * = — -=.GpC F va% \ dx\dx 2 dx 3 \ b\dbib 2 db 2 (j)B{xi,bi) 
v 3 io Jo 

■ {{(1 + x 2 ) <\>k{x 2 )<\>k* (x 3 ) 

+r K (1 - 2x 2 ) [<Pk{x 2 ) - (j) T K {x 2 )] (p K *(x 3 ) ~ 2r K ,(j)j <: (x 2 )(t) s K ,{x 3 ) 

-2r K r K , [(2 + x 2 ) 0£(x 2 ) + x 2 <f K {x 2 )\ (f) s K *{x 3 )} 

■E^{t q ,l 2 )h e {x 2 ,x u b 2 M) 
+ {2r K (p^(x 2 )(t)K*(x 3 ) - 4r^r^.0^(x 2 )0x*(^3)} 

•^e;, 0^1^2,61,62)}, (54) 

where tk = tuq /raB,TK* = TTiK*/fn Bl the evolution factors take the form of 

E^\t,l 2 ) = C^(t,l 2 ) a 2 (t)-exp[-S ab ] } (55) 

with the Sudakov factor S a t, and the hard function h e (xi,x 2 , 61, b 2 ) as given in Eq. (1A9I) 
and Eq. (1A1I) respectively, and finally the hard scales and the gluon invariant masses are 



t q = max( v /i 2 "m jB , y /x 1 x 2 m B , a/(1 - x 2 )x 3 m B , I/61, l/6 2 ); , 

t' q = max( v /x7m jB , y /x 1 x 2 m B , a/|x 3 - Xi\m B , I/61, l/6 2 ), (56) 

/ 2 = (1 - x 2 )x 3 m 2 B - |k 2T - k 3T | 2 « (1 - x 2 )x 3 m 2 B , 

I' 2 = (x 3 - Xi)m 2 B - |kiT - k 3T | 2 « (a? 3 - Zi)to|- (57) 

Finally, the total "quark-loop" contribution to the considered B — » if if* decays can 
be written as 



<. = < if*if l^ls >= £ \ m# jc + m<&. , (5 



r(«) 1 «>f(9) 



where A g = V^,V^. 

From the functions C( q >((j,,l 2 ), one can see that the quark-loop amplitudes depend on 
both the renormalization scale /i and the gluon invariant mass I 2 . In the naive factorization 
approach, the assumption of a constant I 2 , I 2 ~ rn 2 /2, introduces a large theoretical 
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FIG. 6: Chromo-magnetic penguin (Osg) diagrams contributing to B° 
B + _+ x* + K° and K + K*° decays. 



K*°K° + K°K*°, 



uncertainty as making predictions. In the pQCD approach, however, I 2 is related to the 
parton momenta unambiguously. Because of the absence of the end-point singularities 
associated with I 2 , 1' — > 0, in Fig. (j5Ja and fl5}b respectively, we have dropped the parton 
transverse momenta kx in Z 2 , Z' 2 for simplicity [131 ] . 

From Fig. ([5]), it is easy to see that the "quark- loop" diagrams contribute only to 
B° -»> iT* Z° + f^F* and 5+ -> fC*+¥°, if+if decays. For B° -+ K + K*~ and 
K* + K~ decays, there is no such kind of NLO contributions. 



C. Chromo-magnetic penguin contributions 



As illustrated by Fig. (jHJ), the chromo-magnetic penguin operator Og g also contribute to 
B — > KK* decays at NLO level. The corresponding weak effective Hamiltonian contains 
the b — > dg transition, 



nicmp 
n eff 



■^Wd cti f o 8 



with the chromo-magnetic penguin operator, 



O, 



9s 



57T" 



^d^il+^T-G-b 



i] ia> 3 ' 



(59) 



(60) 



where i,j being the color indices of quarks. The corresponding effective Wilson coefficient 
r*ff 



+ c 5 
(ill 
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c * 

In Ref. [22|, the authors calculated the chromo-magnetic penguin contributions to 
B —>■ <pK decays using the pQCD approach. They considered nine chromo-magnetic 
penguin diagrams corresponding to the non-local operator 0' s , as given in Eq. (2.3) of 
Ref. [22J, generated by operator 0% g as defined in Eq. fl60|) . The first two Feynman 
diagrams (a) and (b) in Ref. [22j are the same as Figs. |S|a and (jn])b here. According 
to Ref. [22|, the diagrams (a) and (b) dominate, while other seven diagrams are small or 
negligible. It is therefore reasonable for us to consider the NLO contributions induced by 
the diagrams (a) and (b) only, for the sake of simplicity. 
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The decay amplitude for Figs. Eh and Eb can be written as 

4 f 1 r°° 

Mk*k = —/=GfCp7ti% / dx\dx 2 dx 3 \ b\db^b 2 db 2 (^ B {x\,b\) 
v3 Jo Jo 

■ {{- (1 - x 2 ) [2(f)K*{x2) - r K * [3<f>' K .(x 2 ) - #f.(x 2 )] 
-r K ,x 2 [<f>' K .{x 2 ) + 0a-.(« 2 )]] </4fe) 
+r x (l + x 2 )x 3 -0x*(x 2 ) [30^(x 3 ) + 0^(x 3 )] 
-r K «r K (1 - x 2 ) [0x*(^a) + 0k*(^)] [30^(x 3 ) - 0^(x 3 )] 
-r,K,r x X3 (1 - 2x 2 ) [<^*(x 2 ) - 0x*(z 2 )] [30^(x 3 ) + $c(z s )] } 

■E g {t q )h g {A,B,CMMM^ 2 ) 
+ {4r K *<p s K *{x 2 )<j)£{x 3 ) - 2r K *r K x 3 <j) s K *{x 2 ) [30£(x 3 ) + 0£(z 3 )] } 

■E g (t' q )h g (A' ,B> ',C",Mi,&3,3i)} , (61) 

for the case of i? — > i^* transition, and 

Mkk* = —/=GFCpm% / dx\dx 2 dx 3 \ b\db]b 2 db 2 $ B {x\,b\) 
v3 Jo io 

{{-(1 -x 2 ) [20£(x 2 ) +r K [30£(x 2 ) -</£(x 2 )] 
+r x x 2 [</>£(x 2 ) + 4> T K {x 2 )]] (p K *(x 3 ) 
+r K * (1 + x 2 ) x 3 <f)£{x 2 ) [&f> 8 K .(x 3 ) + <^*(x 3 )] 
+r x r x * (1 - x 2 ) [#£(z 2 ) + 0ife)] [30^»(x 3 ) - ^.(x 3 )] 
+r K r K *x 3 (1 - 2x 2 ) [</£(x 2 ) - $fe)] [30^(s 3 ) + ^*(x 3 )] } 

• J B fl (g/i g (A, J B ) C,fo 1 ,6 2 ,6 3 ,x 2 ) 
- (4r^0^(x 2 )(/)x*(x 3 ) - 2r^r^.x 3 0^(x 2 ) [30^, (x 3 ) + </4»(x 3 )] } 

■^(g/^A', 5', C", 6a, 6i, 6 3 , xi) } , (62) 

for the case of B — > K transition. Here the hard scale t q and t' are the same as in Eq. ( 1561) . 
The evolution factor E g (t) in Eqs. ( IBTi) and ( 162]) is of the form 

£,(*) = C|f(0 «2(*) • exp [-5,] , (63) 

with the Sudakov factor S g and the hard function h g , 

S mg (t) = s (x 1 m B /V2, 61 J + s ( x 2 m B /V2, b 2 ) + s (( 1 - x 2 )m B /V2, b 2 
+s (x 3 m B /V2, b 3 ) +s((l- x 3 )m B /V2, 



1 



# 



ln ln(t/A) + ^ ln(t/A) + ^ ln(t/A) 



ln(6iA) -ln(6 2 A) - ln(6 3 A) 



(64) 



h g {A, B, C, b x , b 2 , 63, x,) = -S t (x 4 ) iTo(S6i) K (C6 3 ) 

/•t/2 

•/ d6tcm6 ■ J {Ab! tern d)J (Ab 2 tan 6) J {Ab 3 tan 6), (65) 

io 



18 



where the functions Kq(x) and Jq(x) are the Bessel functions, the form factor St(xi) with 
i — 1,2 has been given in Eq. (jA7p . and the invariant masses A^\B^ and C^ of the 
virtual quarks and gluons are of the form 

A = v^rriB, B — B' — y/xTx^m B , C = za/(1 - X2)x 3 m B , 
A' = y/xirriB, C = y/x\ — x 3 m B . (66) 

The total "chromo-magnetic penguin" contribution to the considered B —> KK* decays 
can therefore be written as 



M%$> = < K*K\H c ™j\B >= A t 



M (g) , M (g) 

1V1 K*K + 1V1 KK* 



(67) 



where \ t = V th V* d . 

From Fig. ([61) , one can see that the chromo-magnetic penguins contribute only to 
B° -► iT* Z° + K°K*° and 5+ -> K*+K°, K+K*° decays. For 5° -► iT+iT*- and 
K* + K~ decays, there is again no such kind of NLO contributions. 



V. NUMERICAL RESULTS AND DISCUSSIONS 

A. Input parameters 

Besides those specified in the text, the following input parameters will also be used in 
the numerical calculations: 

m B = 5.28GeV, m K = 0.49GeV, m K * = 0.892GeV, m b = 4.8GeV, 

m 0K = 1.7GeV, m w = 80.41GeV, m t = 168GeV, a em = 1/128, 

f B = 0.21GeV, f K . = 0.217GeV, f K , = f K = 0.16GeV, 

t b o = 1.528ps, t b + = 1.643ps, (68) 

For the CKM quark-mixing matrix, we use the Wolfenstein parametrization as given 



in Ref. 23 



V ud = 0.9745, V us = A = 0.2200, \V ub \ -- = 4.31 x 10" 



/, | V U0 | ^.KJA. XN i^ , 



K, = -0.224, V cd = 0.996, V cb = 0.0413, 
|Kd| = 7.4 x 10- 3 , \/ 4s = -0.042, \V tb \ = 0.9991, (69) 

with the CKM angles f3 = 21.6°, 7 = 60° ± 20° and a = 100° ± 20°. The unitarity 
condition V ub V* q + V cb V* q + V tb V t * q = for q = d, s is employed 

B. Branching ratios 

In the pQCD approach, the form factor A^ K (q 2 = 0) and F^{* K (q 2 = 0) can be 
extracted from the decay amplitude F e x* and F e x* as shown in Eqs. (fl9~|) and fl30l) . via 
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the following relations, 



^B^Kf „2 _ n \ _ * eis: 



F " (9=0) = sS' (70) 

^* V ^W (71) 

Consequently, one can find the NLO pQCD predictions for the values of the corresponding 
form factors at zero momentum transfer: 

A B-+K*^ q 2 = o) = 0.38 ± 0.05(w 6 ), Fg{* K (q 2 = °) = °- 36 ± 0.06(cj 6 ), (72) 

for uif, = 0.40 ± 0.04GeV, which agree well with those obtained in QCD sum rule calcula- 



tions, for example, in Refs. [18] , ll 



For a general charmless two-body decays B —* f with / = M 2 M 3 , the branching ratio 
can be written in general as 

WB - « = T * vL; }Mf (73) 

where tb is the lifetime of the B meson, and M. =< KK*\TC e ff\B > for the case of 
f = KK*. 

Using the wave functions and the input parameters as specified in previous sections, 
it is straightforward to calculate the branching ratios for the considered decays. For 
B + — > K + K and B + — > K* + K decays, we show in Table HT1 the CP-averaged branching 
ratios 

Br(B ->/) = £ pr(B -> /) + fir(5 -> J)] . (74) 

For S° decays, it is a little complicate since both B° and B can decay into the final 
state / and / simultaneously. In Table HT] we show the CP-averaged Br's for B° —>■ fi, 
B° —> fi and for B° —>■ fi+f\ with f\ = K°K , respectively. The third result corresponds 
to the measured upper limit. For B° — > / 2 , f'2 and B° — > f 2 + fa with / 2 = K + K*~ , we 
take the same convention. 

Except for the LO results, we always use the NLO Wilson coefficients in the calcu- 
lations. The label +VC, +QL, +MP and NLO denote the pQCD predictions with the 
inclusion of the vertex corrections only, the quark loops only, the magnetic-penguin only, 
and all the considered NLO corrections, respectively. For the sake of comparison, we also 
show currently available experimental results [24j and the numerical results evaluated in 
the framework of the QCD factorization (QCDF) Q. 

It is worth stressing that the theoretical predictions in the pQCD approach have rel- 
atively large theoretical errors induced by the still large uncertainties of many input 
parameters. The pQCD predictions for the branching ratios with the consideration of 
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TABLE II: The pQCD predictions for the branching ratios (in unit of 10™ 7 ). The label LO means 
the leading-order results, and +VC, +QL, +MP, NLO mean the predictions with the inclusion 
of the vertex corrections, the quark loops, the magnetic-penguin, and all the considered NLO 
corrections, respectively. 



Mode 


LO 


+VC +QL +MP NLO 


Data 


QCDF 


B+ - 
B+ - 


► K*+K° 


4.2 
2.0 


5.3 

2.7 


5.8 
2.3 


3.1 
1.6 


3.2 
2.1 


< 11 


3.0+|° 
3.0+™ 


B ^ 
B ^ 
B ^ 


K°K* U 
K°K*° 
K°K*° + K°K*° 


2.1 
6.4 
13.7 


3.0 
6.9 

14.0 


2.9 

8.0 
15.2 


1.8 
4.3 

6.7 


2.4 
4.9 
8.5 


< 19 


2.6+J;* 
2.9+™ 


B ^ 
B ^ 
B ^ 


K+K*' 

— *+ 
K-K 

K+K*~ + K~K* + 


1.1 

0.41 
2.7 


— 


— 


— 


0.83 

0.17 

1.3 




n -,4+1.07 

n -,4+1.07 
u - ±4 -0.14 



major uncertainties are the following (in unit of 10 

T*0 



Br{ B+ - 

Br( B + - 

Br( B° 

Br( B° 

Br(B° 

Br( B° - 

Br( B° - 

Br(B° 



-K+K 

■ K* + K° 

♦ K°T° 

>K°K*° 

*fi + fi 
K + K*- 

K-K*+ 

>f2 + f2 



3 -2_ ; 6 (^b) - O '.i(oi) -o'A^k) -0.l( a iK*), 

2A±0.2(uJ b 



-0.0/ 
-0.1 



a 



-0.3/ n 



0S (n 
0AK a iK* 



4.9±S3(^)±83(«)±g5(aur)±85(auc*), 

8.5+?; 2 7 (^) ± 0.1(<x)^ZM±ki(<hK* 



0.83 



-0.04 

-0.08 ' 

0.02 

-0.01 



u Jb )±0.Q(a) + _° f 1 (a iK ) + _° f 8 (a iK , 



0.17±°;°iH) ± 0.07(a)+^(a^) ± 0.02(a^ 

1.3±0.1(^)+°;°( 



-0.0/^N+0.2/ \+0.4/ 

I Oi ) _ 06 [CliK ) _ .3 \fliK* 



(75) 

The major theoretical errors are induced by the uncertainties of oOb = 0.4 ± 0.04 GeV, 
a = 100°±20°, and Gegenbauer coefficients a\K — 0.17±0.17, a 2 K = 0.115±0.115; a\K* — 
0.03 ± 0.03, ci2K* — 0.11 ± 0.11, respectively. Additionally, the final-state interactions 
remains unsettled in pQCD, which is non-perturbative but not universal. Fortunately, 
good a gree ment between the pQCD predictions for the branching ratios of B — > KK 
decays 25j and currently available experimental measurements 24j indicates that the 
FSI effects are most possibly not important. 

From the numerical results, it is easy to see that 

• The pQCD predictions for the branching ratios of B+ — > K+K and B° — > 
K K + K°K*° are consistent with currently available experimental upper limits. 
Inclusion of the NLO contributions decreases the central value of the LO predictions 
by about 30% to 80%. The chromo- magnetic penguin provide the dominant NLO 
contributions. 

• For B° — > K+K*~ decay, the pQCD prediction is rather different from that from 
the QCD factorization approach. Such difference could be tested in the forthcoming 
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LHCb experiments. For other decays, the pQCD predictions agree well with the 
corresponding QCDF results within one standard deviation. 



C. CP- violating asymmetries 

Now we turn to the evaluations of the CP- violating asymmetries of B — ► K*K decays 
in pQCD approach. For B + — ► K + K and B + — ► K* + K decays, the direct CP-violating 
asymmetries Aqp can be defined as: 

Afp = S/i'-i^;: , (76) 



I7W/I 2 + \M f 



2 



where .A/f/ =< f\H e ff\B > and .M/ =< f\H e ff\B >. 

The pQCD predictions for the direct CP-violating asymmetries of the considered decays 
are listed in Table IIHI For comparison, we also reproduce verbatim the corresponding 
numerical results evaluated in the framework of the QCD factorization (QCDF) [201 ] . 



TABLE III: The pQCD predictions for the direct CP asymmetries of B — » KK* decays (in 
units of percent). 



Mode 


LO +VC +QL +MP NLO 


QCDF 


A d S r P (B + - 


-> K + K*") 


-53.4 -39.2 -5.1 -49.8 -6.9 


_24+ 28 
Z4 -39 


A d S r P (B + ~ 


-> K*+K u ) 


8.1 -12.3 7.2 10.6 6.5 


io+29 
1J_ 37 



The pQCD predictions for Afjp and the major theoretical errors for B A 
*o 



1 dir 

K+W" , K* + TT decays are 

ldir(T>+ _^ T^+T"*°\ - r R n+5-6/, , N+l.O/' „,\+9.2/„ .. _\+4.0/ ... \1 v i n -2 



^p(^ + - ^ + = L-6-9i^(o; 6 )^(«)^(a^)^o(a^)J x 10^, (77) 

^(B+ - ir* + Z°) = [6.5+^(a; 6 )+li(a)+^(a^)+ 2 3i(a^)] * 10" 2 , (78) 

where the dominant errors come from the variations of u;& = 0.4±0.04 GeV, a = 100°±20°, 
and Gegenbauer coefficients a\K = 0.17 ± 0.17, a^K = 0.115 ± 0.115; clik* = 0.03 ± 0.03, 
a 2 K* — 0.11 ± 0.11, respectively. 

In Fig. [7J we show the a— dependence of the direct CP-violating asymmetries Afjp for 

B + — > K + K (the solid curve) and B + — > K* + K (the dotted curve) decay, respectively. 
The left figure is for the LO pQCD predictions and the right one for the NLO pQCD 
predictions. One can see from the numbers and figures that (a) as usual, there exist 
relatively large differences between the pQCD and QCDF predictions; (b) the LO and 
NLO pQCD predictions for the direct CP-violating asymmetries are also rather different; 
and (c) the NLO contribution from the "Quark-loops" ("Vertex corrections ") leads to 
the dominate change of Afj P for B+ -► K+K*° ( B + -* K*+K°) decay. 

We now study the CP-violating asymmetries for B°/B — > K + K*~ (K~ K* + ) decays. 
Since both B° and B can decay into the final state K + K*~ and K~K* + , the four time- 
dependent decay widths for B°(t) -> K+K*~, B°(t) -► R-R* + , B°(t) -> K~K*+ and 
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FIG. 7: The direct CP asymmetries (in percentage) of B + — ► K + K (solid curve) and B + — > 
K* + K (dotted curve) as a function of CKM angle a. (a) shows the LO results, while (b) shows 
the NLO results. 



B it) — > K + K* can be expressed by four basic matrix elements: 

g = (K + K*-\H eff \B°), h = (K + K*-\H eff \B°). 
g = {K-K* + \H eff \~B*), h = (K-K* + \H eff \B°) 



(79) 



-rrO 



which determine the decay matrix elements of B° — > K+K* , B — ► K K* + , B° 



-rrO 



K K* + and B — > K+K* at t — 0. Besides the matrix elements g,g, h and h, one also 



-rrO 



need to know the CP-violating parameter coming from the B° — B mixing: 

B 1 =p\B°)+q0 i ), B 2 =p\B°)-q\B°), (80) 

with \p\ 2 + \q\ 2 = 1. 



Following the notation of Ref. [10], the four time-dependent widths are given by the 
following formulae: 

T(B°(t) -> K+K*-) = e- Tt -{\g\ 2 + \h\ 2 ) x {1 + a e > cos(Amt) + a e+e , sin(Amt)} , 
T(B°(t) -> K~K*+) = e- r *i(|7i| 2 + \h\ 2 ) x {1 - at cos(Ami) - a e+ - e > sin(Amt)} , 



T{B°(t) -> K~K*+) = e- rt -{\g\ 2 + \h\ 2 ) x {1 - a? cos(Amt) - a t+ ? sin(Amt)} , 
r(S°(*) -> K+K*-) = e- Tt h\g\ 2 + \h\ 2 ) x {1 + a e > cos(Amt) + a e+e , sin(Amt)} ,(81) 



where the CP -violating parameters are 

\g\ 2 -\h\ 2 

\g\* + \h? 

2 -\g\ 2 



a r / 



o 7 ' 



\M 

\h\ 2 + \g\ 



flf-i-?' 



-2/m(J|) 

i + IW 

v p h> 

1 + \g/h\ 2 



(82) 
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FIG. 8: The CP-violating parameters of B°/B° -> K K*°(K°K* ) decays, (a); and B°/B° -► 
K + K*~ (K~ K* + ) decays, (b) : a e i (dash-dotted line), a^ (dotted line), a e+e / (dashed line) and 
a-e+T' (solid line) as a function of CKM angle a. 



with q/p = e 2i/3 a and /3 = 21.6° is one of the three CKM angles. 
Similarly, the four time-dependent decay widths for B° — > 



K Q K*°, B° 



B° — > K K*° and B — > K K*° (here K° means Kg) can also be defined as 



-*o, 



-*o, 



,-rrO, 



(K°K*"\H eff \B°), h = (K°K™\H eff \W), 



g = (irV |^ e// |5°), h = (K"K w \H eff \B u ), 



K°K 



■*o 



33) 



One can define, consequently, the four CP-violating parameters <v, a e+e >, a? and a e+I > for 
B°/B° — > fi + fi decays in the same way as in Eq. (!82|) . In Fig. [HJ we show the pQCD 
predictions for the eight CP-violating parameters for the considered decays. 

The central values of the pQCD predictions for the CP-violating parameters are 



a e i 



0.13, 



-0.96, 



<v 



-0.72, 



0.59. 



54) 



for B°/R -> K+K*- + K~ K* + decays, and 



a, = -0.50, 



K°K*° + K°K*° 



0.24, 



a T i 



0.05, 



0.12, 



15) 



for B°/W -* K U K" + K"K* U decays. 

As pointed in Ref. 26j, it may be conceptually incorrect to evaluate the Wilson coeffi- 
cients at scales down to 0.5 GeV. The explicit numerical values for the Wilson coefficients 
Ci([/,) — Ci (/i) for /j, = 0.5, 1.0, 1.5 and 2.0 GeV, as listed in Table HVj also support this 
expectation: the values of the Wilson coefficients Cs^^fi) at // = 0.5 GeV are about 
four to seven times larger than those at n — 1.0 GeV. For C 5 (fi), specifically, C 5 (0.5) 
and C5(1.0) even have a different sign besides the large difference in their magnitude. In 
the region of /i > 1.0 GeV, however, the //—dependence of all Wilson coefficients become 
relatively weak. It is therefore reasonable for us to choose //q = 1-0 GeV as the lower 
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cut-off of the hard scale, instead of /io = 0.5 GeV as being assumed in Ref. [13j. We then 
fix the values Cj(/i) at Cj(/io = 1-0), whenever the scale \i runs to below the scale /i - 

In order to show directly the /zo-dependence of the branching ratios and CP- violating 
asymmetries, we recalculated these quantities for B — > KK* decays by setting /i = 0.5, 
1.0, 1.5 and 2.0 GeV, respectively. It is easy to see from the numerical results as listed in 
Table |V] that the pQCD predictions are relatively stable against the variation of /io for 
/i > 1-0 GeV. We therefore set /io = 1-0 GeV to be the cut-off scale for Wilson coefficients 
Ci(fj,). Of course, the issue of /i -dependence need more studies. 



TABLE IV: NLO Wilson coefficients C;(/z) for /i = 0.5 - 2.0 GeV, respectively. 



Mo 


Ci(n) C 2 (m) 


CM 


CM) 


C 5 (H) 


CM 


CM CM 


cm dM 


0.5 GeV 


-0.9923 1.6537 


0.1729 


-0.3122 


-0.1143 


-0.8276 


0.0010 0.0056 


-0.0148 0.0092 


1.0 GeV 


-0.5093 1.2790 


0.0428 


-0.0898 


0.0150 


-0.1321 


-0.0002 0.0008 


-0.0120 0.0050 


1.5 GeV 


-0.3773 1.1920 


0.0289 


-0.0652 


0.0153 


-0.0856 


-0.0002 0.0005 


-0.0112 0.0038 


2.0 GeV 


-0.3114 1.1518 


0.0230 


-0.0541 


0.0145 


-0.0672 


-0.0002 0.0004 


-0.0108 0.0032 



TABLE V: The pQCD predictions for the branching ratios (in unit of 10 7 ) and direct CP- 
violating asymmetries (in unit of IO -2 ) for the considered B — ► KK* decays, assuming /io = 
0.5, 1.0, 1.5 and 2.0 GeV, respectively. 



Mode 


/x = 0.5 


/io = 1.0 


Ho = 1.5 


/i = 2.0 


Br{B+ -» K + K*") 
Br(B+ -► K*+K°) 
Br(B°/B° - A + A) 
Br(B°/B° -> h + h) 


4.7 
2.6 
22.5 
5.4 


3.2 

2.1 
8.5 
1.3 


2.6 
1.3 

5.0 
0.78 


2.1 

0.8 
3.5 
0.55 


A d S r P {B + -» K + K*") 
A% r P (B + -» K*+K°) 


-4.0 

16.8 


-6.9 
6.5 


-7.1 
-1.5 


-5.1 
-5.8 



VI. SUMMARY 

In this paper, we calculate some NLO contributions to the branching ratios and CP- 
violating asymmetries of B — > KK* decays in the pQCD factorization approach. 

From our calculations and phenomenological analysis, we found the following results: 

• The NLO contributions from the QCD vertex corrections, the quark-loops and the 
chromo-magnetic penguins can be rather large and provide significant modifications 
to the LO predictions. 
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The NLO pQCD predictions for the form factors of B — ■*■ K* and K transitions are 

A*-* K \(f = 0) = 0.38±0.05(cj 6 ), 
F o*,r K (Q 2 = 0) = 0.36 ± 0.06(w 6 ), (86) 

for Ub = 0.40 ± 0.04GeV, which agree well with those obtained in QCD sum rule 
calculations. 

The pQCD predictions for the branching ratios are 

Br {B+ -► K*°K + ) = 3.21JJ x lO" 7 , 
5r(5+ -► K* + K°) = 2.1±\i x 10~ 7 , 



Br(B° -► K K*° + K°K* ) = 8.511;? * lO" 7 , 
Sr(S° -»• A'+A'* + K-K* + ) = 1.3lg;| x 10~ 7 , (87) 

where the theoretical errors from various sources are added in quadrature. These 
pQCD predictions are consistent with both the QCDF predictions and currently 
available experimental upper limits. 

The direct CP- violating asymmetries for B + — ► K + K , K* + K are (in unit of 10~ 2 ) 

.4&(B + - K+iT ) = -6.91J5;|, 



^(5+ -► 7T+7T) = 6.5™, (88) 

which are rather different from those in the QCDF approach. 
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APPENDIX A: RELATED FUNCTIONS 

We show here the function h^s, coming from the Fourier transformations of H^°\ 



h e (x 1 ,x 2 , h, b 2 ) = K ( v /x 1 x 2 m B b 1 ) [$(bi - b 2 )K (^/x^m B bi) I {^fx^m B b 2 ) 

+9(b 2 - b^Ko {^/x^m B b 2 ) I (y/x^m B bi)] S t (x 2 ), (Al) 



h a (x 2 , x 3 , b 2 , b 3 ) = K (i\/(l - x 2 )x 3 m B b 2 ) [6(b 3 - b 2 )K (iy/x^mnbz) I {i^/x^m B b 2 ) 

+6(b 2 - b 3 )K {i^/x^m B b 2 ) I [i^/ximsb^] S t (x 3 ), (A2) 
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h f (xi,x 2 ,x 3 ,b 1 ,b 3 ) = < 9(bi - h)K (mBy / xiX2bi)lo(mBy / x 1 x 2 b 3 ) 

+ 6{b 3 - b 1 )K (mBv / x 1 X2b 3 )I (m B y/x 1 x 2 b 1 ) 

f L H ( v /(x 2 (x 3 - xi))m B b 3 ), for x 1 - x 3 < . , , , 
K^ytt^m^a), for x 1 -x 3 >0'- (Ao) 



^(xi,X2,X3,6i,6 3 ) = < 6>(6i - 6 3 )K («v / ( 1 ~ zj^&iM m )I (z\/(l - x 2 )x 3 b 3 m B ) 



+ (0(6 3 - 6i)K (i\/(i - x 2 ) x 3 637715)10(2 a/ (! - x 2 )x 3 bim B ) 

K (m By /{ 3c 1 - x 3 ){\ -~a^)6i ), for xi - x 3 > \ .^^ 
f H^KV^ - ii)(l - s 2 )&i), for X! - x 3 < ) A 



/i}(xi,X2,x 3 ,6i,6 2 ) = < 6>(6i - 6 3 )K («v / ( 1 - x 2 )x 3 6im j B)Io(iA/(l - x 2 )x 3 b 3 m B ) 

+ #(6 3 - 6i)K (z-/(i - x 2 )x 3 6 3 M i j)Io(V( 1 - ^2)x 3 6im B ) 

K (m B Fi6i), for F? > ()\ , , 

-^\m B ^\Ff\bi), for F x 2 < Oj ' l J 

where Jo is the Bessel function and X0, -^0 are rnodified Bessel functions K Q (—ix) = 
— (tc/2)Y (x) + i(ir/2)Jo(x), and F^s are defined by 

if 1} = l-^a-Xa-xO. (A6) 



The threshold resummation form factor St(xi) is adopted from Ref.[27j. It has been 
parametrized as 

* w = 2 ^y mi -*)]•, <A7) 

where the parameter c = 0.3. 

The evolution factors Ee and M appeared in Eqs. ( fl9l) to ( 1381) are given by 



£ e (£) = <*,(*) exp^a^i)], 

E' e (t) = a s {t)eKjp[-S cd {t)]\ b2=bl , 

E a {t) = a s {t)exp[-S gh (t)}, 

Kit) = a s (t)exp[-S ef (t)]\b 2 =b 3 , (A8) 
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where the Sudakov factors can be written as 

S a b(t) = s [xim B /V2, foij + s (x 2 m B /V2, b 2 ) + s Ml - x 2 )m B /V2, k 
1 



In Ht / A l I In Ht/A) 



SM 



-ln(6xA) -ln(fo 2 A) 

s (xim B /V2, b\ J + s (x 2 m B /V2, bA + s ( (1 

+s (x 3 m B /V2 1 b 3 ) +sf(l- x 3 )m B /V2, b 3 
1 



x 2 )m B 



/V2,h 



Pi 



21nJ« + ln>W A > 



ln(foiA) ' ~-ln(6sA). 
S e /(£) = s \xi171b/V2, bij+s (x 2 m B /V2, b 2 )+s((l- x 2 )m B /V2, b 2 j 

+s (x 3 m B /V2, b 2 ) + s (( 1 - x 3 )m B /V2 : M 

_J_r ln(t/A) ln(t/A) - 

A L -lB(6iA) +zm -lB(62A). ' 

^(i) = s ( x 2 m B /V2, b 2 \ + s (x 3 m B /V2, b 3 ) +s((l- x 2 )m B /V2, k 
+ s 



In ln(tM) , I In ln(t/A) 



-ln(6iA) -ln(6 2 A) 



X Z )m B />/2,kij - -5- 



where the function s(g, 6) are defined in the Appendix A of Ref. [6]. 
The hard scale ij's appeared in Eqs. f|T9|) to fl38l) are of the form 



£„ = max( v /x 2 "m jB , ^/x^tub, I/61, l/6 2 ) , 

t' a = max( v /x7m jB , y /x 1 x 2 m B , 1/6 1( l/6 2 ) , 

t b = max( v / x 2 |l - x 3 - ^Jm B , y/x 1 x 2 m B , l/&i, I/63) , 

4 = max(A/x 2 |x 3 - Xi|m B , ^XiX 2 m B , l/&i, I/63) , 

t c = max(v(l - x 2 )x 3 m B , y/\x[ - x 3 \(l - x 2 )m B , I/61, 1/63) , 

t(. = max( v / |l - 2:2(1 - £3 - xi)|m s , a/(1 - x 2 )x 3 m B , l/b u l/b 3 ) 

t d = max( v / (l -x 2 )x 3 m B , a/(1 - x 2 )m B , l/b 2 , 1/63) , 



max 



(a/(1 - x 2 )x 3 m B , y/x^m B , l/b 2 , 1/h). 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 
(A14) 
(A15) 
(A16) 
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